Flux lattice melting and depinning in the weakly frustrated 2D XY model 



S0ren A. Hattel and J. M. Wheatley* 
Interdisciplinary Research Centre in Superconductivity, University of Cambridge, 
Madingley Road Cambridge, CBS OHE, United Kingdom. 
( February 1, 2008) 



in 

On 
OS 

PL, 



> 

in 
in 
o 

<N 
O 
in 
o> 



C 

O 
o 



Monte Carlo simulations of the frustrated 2D XY model were carried out at small commensurate 
values of the frustration /. For / = 1/30 a single transition was observed at which phase coherence 
(finite helicity modulus) and vortex lattice orientational order vanish together. For / = 1/56 a 
new phase in which phase coherence is absent but orientational order persists was observed. Where 
comparison is possible, the results are in detailed agreement with the behavior of the lattice Coulomb 
gas model of vortices. It is argued that the helicity modulus of the frustrated 2D XY model vanishes 
for any finite temperature in the limit of weak frustration /. 
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Motivated by high temperature superconductors, there 
has been renewed interest recently in the nature of phase 
transitions, if any, in strongly itype II 2D superconduc- 
tors in applied magnetic neldsa. One approach to the 
2D problem in the infinite k limit is provided by the 
lowest Landau level approximation to Ginzburg-Landau 
theory, where a truncated basis consisting only of states 
in the lowest Landau level of linearized Ginzburg-Landau 
theory is retainedo. First order melting transitions have 
been observed in Monte Carlo simulations using this 
approachla, although_no transition was observed in sim- 
ulations on a sphereu. An alternative discretization pro- 
cedure places Ginzburg-Landau functional on a lattice, 
usually by use of phase only models such as frustrated 
XY or Villain models. Here fluctuations in the amplitude 
of the superconducting order parameter are neglected, 
which is a reasonable approximation far from the critical 
temperature. In the XY and Villain models compactness 
(equivalence of states whose local superfluid phases differ 
by multiples of 2tt) has been imposed in different ways. 

The interaction potential of the Villain model has very 
special features. The model shows a separation of vortex 
and spin wave degrees of freedom; namely, after dual- 
ity transformation, the partition sum factorizes into a 
product of spin wave (gaussiaa) fluctuations and 2D lat- 
tice Coulomb gas components! The lattice Coulomb gas 
amounts to a 2D "lattice London model" for the vortex 
state. It was on the basis of the continuum Coulomb 
gas model that Huberman and DoniachQ and FishenJ dis- 
cussed the vortex lattice melting transition in 2D super- 
conductors. Here Kosterlitz-Thouless dislocation medi- 
ated melting theoryma leads to a melting temperature es- 
timate which is independcnt-pf vortex density. However, 
in a series of papers Moorell3, has implicitly questioned 
the validity of the London model, and argued that the 
lower critical dimension for superconductivity is in fact 
d = 4. 

The question naturally arises whether the phase only 



discretization of Ginzburg-Landau theory provided by 
the Villain model is fully equivalent to other discretiza- 
tion schemes. If it were not, then the validity of conclu- 
sions based on the Coulomb gas (London) model would 
be in doubt. For example, in the alternative discretiza- 
tion provided by the XY model, the separation of vortex 
and spin- wave degrees of freedom is not complete. Thus 
the London model is not strictly valid even for the vor- 
tex contribution to the thermodynamics. Both the lattice 
London and XY models have been extensivelyO used in 
studies of vortex lattice melting, and it is essential to 
know if they indeecLshare the same continuum limit. 

In a recent LetterEa Franz and Teitel reported very de- 
tailed Monte Carlo simulations of the lattice Coulomb gas 
(LCG) in two dimensions. While the discreteness of the 
lattice is expected to become unimportant in the dilute 
limit, these authors found that the "continuum limit" is 
subtle, and only reached for surprisingly small values of 
the density (denoted /). In sufficiently dilute systems, a 
"depinning" or "floating" phase transition occurs from a 
low temperature phase with long range translational or- 
der (LRTO) to a phase with algebraic order, followed by 
a first order 2D melting transition to a disordered phase. 
For dense systems, the phase with algebraic order is ab- 
sent, i.e. depinning and melting transitions coincide. The 
existence of a phase with LRTO is an artifact of pinning 
in the lattice model because there is a finite energy cost 
for displacing a vortex in the ground state configuration. 
The appearance of a phase with algebraic order signals 
the onset of the continuum limit. 

The transition temperature T p at which LRTO disap- 
pears vanishes linearly in /. The 2D melting tempera- 
ture, where algebraic order disappears, is independent of 
density as anticipated from 2D melting theory The lin- 
ear dependence of T p on density ip-.the dilute Coulomb 
gas can be established analyticallyEl The 2D harmonic 
Coulomb solid on a lattice can be mapped to a dislo- 
cation (vector plasma) problem, which has a continuous 
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FIG. 1. The configuration of vortices in the ground states 
for (a) / = 1/30, and (b) / = 1/56. For convenience only 
system sizes 24 x 24 are shown. 



Kosterlitz-Thouless-Halperin-Nclsorolj unbinding transi- 
tion, provided that the core energy for a dislocation is 
sufficiently large. In the dilute limit the core energy 
can be shown to grow as log(l//) so that this condition 
is always satisfied. In the low temperature phase dis- 
placement fluctuations of the Coulomb solid are screened 
which permits LRTO even in 2D. 

In this paper we report the results of a Monte Carlo 
study of the 2D XY model in the limit of weak commen- 
surate frustration (/ = q = integer) with periodic 
boundary conditions. To avoid biasing the vortex system 
towards hexagonal order, we choose a square background 
grid. The 2D XY model is: 
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where 9i is the phase of the superconducting order pa- 
rameter defined on an L x L square lattice of points la- 
beled by i. The link field Aij (line integral of the vec- 
tor potential between nearest neighbor lattice sites i and 
j) satisfies JZ-Ay = 27r<I>/$o = 2irf. The sum is taken 
counter clockwise around a plaquette. The magnetic flux 
per plaquette $ is uniform which effectively imposes the 
infinite n limit. $0 is the flux quantum, and / is the frus- 
tration. The local vorticity v(R) takes values 0, ±1 and 
is defined via mod(#i — 8j — AjJ = 2w(f—u) where the 
sum is around the plaquette Ro. mod(ir) adds integer 
multiples of 2-7T to x to bring it into the range (— vr,7r]. 
With periodic or anti-periodic boundary conditions the 
total vorticity is constrained to equal the total applied 
flux; i.e the plaquette sum J^rK-^-) = L 2 f . 

The ground state spin configurations of the XY model 
in rational (/ = p/q) magnetic fields change in a highly 
non-trivial way with j>|-and q and have been the sub- 
ject of a large literatureE.3. We have succeeded in finding 
ground states for / = 1/30 and / = 1/56 appropriate to 
this study by slow cooling of q x q unit cells, using peri- 
odic boundary conditions for / = 1/30 and mixed peri- 
odic and anti-periodic conditions for / = 1/56. The vor- 
tex lattice is nearly triangular in these ground states as 
shown in Fig. [I]. However, strictly speaking these states 
have only two-fold symmetry. Near hexagonal symmetry 
is associated with extra stability of these ground states. 

At finite temperatures we measure the helicity modu- 



lus and six-fold orientational order parameter. The helie=. 
ity modulus, T, is equivalent to the superfluid density££l 
and is a measure of long range phase coherence. It is 
defined as the sensitivity of the free energy to a twist in 
the boundary condition along a particular direction via 
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where <5 is the twist ang lei. Note that, for periodic 
boundary conditions, T measures the free energy shift 
in the presence of a flux loop of strength S about the in- 
terior of the torus. The helicity modulus is therefore a 
gauge invariant response. 

In the continuum, a 2D solid has finite orientational 
order below it's melting temperature. The six-fold orien- 
tational order parameter is, 
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where <fik is the angle between a fixed direction in the 
XY plane and the direction of the bond between vortex 
k and its nearest neighbor. tp$ = 1 for a perfect triangular 
lattice. 

We used a heat bath method for the simulations heat- 
ing up slowly from the ground states. For / = 1/30 we 
simulated systems with linear dimension L = 60and90, 
and for / = 1/56 we simulated systems of size L = 56 
and 112. Before computing any averages the system 
was carefully thermalized at each temperature discard- 
ing 120,000 - 680,000 Monte Carlo sweeps (MCS) over 
the entire lattice. From these thermalized states we com- 
puted four averages using 20, 000 MCS, from which the fi- 
nal averages were calculated and the error estimated from 
the standard deviation. For the simulations we used sev- 
eral months of CPU time on a HP7000 workstation. The 
task is more difficult than simulations of the LCG, since 
there are 1/ / times as many degrees of freedom involved 
in the updating algorithm. Furthermore, the system sizes 
accessible are restricted by the fact that the .ground states 
are only periodic on at best 5x9 systemso. For these 
reasons careful finite size scaling was impractical. 

A second important difference between XY and LCG 
models from the point of view of Monte Carlo simulation 
is the presence of energy barriers to the motion of isolated 
vortices in the XY model. To move a vortex through 
one lattice constant requires finite motion of spins in the 
vicinity of the vortex. This energy barrier has been es- 
timated to be E p ~ 0.2(13, and plays an important role 
in the physics of Josephson junction arrays. The bar- 
rier means that equilibration times become long at tem- 
peratures T <C E p . Surprisingly, and for reasons not 
completely clear to us, this problem was less severe than 
expected, and the system equilibrates even at tempera- 
tures of order 0.1E p , as evidenced by our ability to find 
ground states by cooling. Further details of the Monte 
Carlo procedure and additional results from these inves- 
tigations will be presented elsewhere^. 
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FIG. 2. Orientational order and helicity modulus for (a) 
/ = 1/30, L = 60 and L = 90, (b) / = 1/56, L = 56 and 
L — 112. Errorbars correspond to ± one standard deviation. 
Lines are guides to the eye. 

Fig. H shows the computed curves for the T and (p§. 
For / = 1/30 orientational order and long range phase 
coherence remain finite up to T p ss 0.045, where they 
both drop rapidly to values around zero. These data are 
consistent with a single T p at which both order parame- 
ters vanish. For / = 1/56 the helicity modulus vanishes 
at T p w 0.03. At this temperature the orientational order 
also falls, but to a clearly finite value. Orientational order 
persists up to T m ks 0.05. The finite size effect is small. 
Note that the melting temperature has not changed sig- 
nificantly between / = 1/30 and / = 1/56. 

As discussed above, to fully characterize the three 
phases and extract exponents for the two phase tran- 
sitions is a difficult task within the framework of XY 
model simulations. However, it is possible to identify 
the phases with reasonable certainty by analogy with the 
LCG, and by qualitative features of the structure func- 
tion. The results for / = 1/56 clearly show the existence 
of two distinct phase transitions. The first one at Tp is 
associated with the depinning of the vortex lattice from 
the underlying mesh. The helicity modulus vanishes at 
the depinning transition, but orientational order persists, 
consistent with the idea that the vortices form a lattice 
with at least algebraic translational order. 

These-, results should be compared with Franz and 
Teitel'al3 LCG simulations on a square lattice, who found 
that T p < T m for / < 1/30 Their / = 1/60 results are 




FIG. 3. Vortex structure function in the square lattice Bril- 
louin zone for / = 1/56 and L = 112 at temperatures (a) 
T = 0.02, (b) T = 0.04 and (c) T = 0.06. 



equivalent, in our units, to T p = 0.028 and T m = 0.047. 
Comparing Fig. |^ we see there is semi-quantitative agree- 
ment between the behavior of the LCG model and the XY 
model. In particular comparison with Fig. 1 of ref. [12] 
shows that there is a striking similarity between the be- 
havior of the inverse dielectric function of the LCG and 
the helicity modulus of XY model. While the connection 
between e _1 (T) and T(T),ean be rigorously established 
only for the Villain modeB, the close analogy between 
these two quantities is evident. 

This interpretation is supported by the behavior of the 
vortex lattice structure function, 



= ;^E ex P[^- R *]MR>(0)}. 



(4) 



R, 



Fig. ^ shows <S(q) for / = 1/56 for temperatures 
T = 0.02, 0.04 and 0.06 corresponding to the three phases 
of Fig. |[b. In the low temperature phase the Bragg peaks 
are well defined out to high orders. In the intermedi- 
ate temperature phase, the Bragg peak structure is still 
present, but the intensity of each peak is much weaker, 
and falls quickly with order. The near hexagonal sym- 
metry remains clear, compatible with algebraic order. In 
the high temperature phase the Bragg peak structure is 
absent and only the circular symmetry characteristic of 
a liquid is present. 

Long range phase coherence disappears at the depin- 
ning transition, leading directly to the conclusion that 
the unpinned ideal vortex lattice with algebraic ordcxJias 
no long range phase coherence and is not superfluidEEltiJ. 
Moreover, on the basis of the close analogy between the 
LCG and frustrated XY models, T p is expected to vanish 
linearly in the limit of small / and thus the superfluid 
stiffness of the XY model vanishes in this limit. This is 
remarkable because, as is well known, the helicity modu- 
lus of the unfrustrated XY model is finite up to the vortex 
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unbinding transitions at T c = 0.9c3. Here the transition 
mechanism is the unbinding of vortex-antivortex pairs 
and the helicity modulus shows a finite universal jump. 
The origin of singular behavior in the small / limit is 
that a depinned Coulomb gas gives rise to metallic screen- 
ing of the 2D Coulomb vortex-anti- vortex interaction at 
long lengthscales, which means that thermally excited 
free vortices appear at any finite temperature, and the 
vortex unbinding transition is absent. 

It would be interesting to know the critical value of the 
commensurate frustration, such that for all / < f c ,T p lies 
below T m . However it appears unlikely to us that T m — 
T p for any / < 1/56 because / = 1/56 is particularly 
strongly pinned, since it has near hexagonal symmetry. 
It is also possible that depinning and melting transitions 
separate for some values of / > 1/30. While from our 
simulations we cannot rule out the presence of an hexatic 
phase for T > T m , this was not observed in the LCGI3. 

In conclusion, the lattice Coulomb gas and frustrated 
XY models show qualitatively very similar melting and 
depinning phase diagrams in low vortex density (weak 
commensurate frustration) limit. In complete agreement 
with the arguments of Mooral3, we find that thermal fluc- 
tuations of the vortex lattice destroy long range phase 
coherence in the continuum limit of the frustrated XY 
model at any finite temperature. On the other hand, 
the model shows a genuine thermodynamic phase tran- 
sition, which can be identified with the disappearance of 
algebraic long range order in the gauge invariant current 
pattern corresponding to the vortex lattice structure. 
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